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ABSTRACT 

The correlation between the cosmological rest-frame v¥ v spectrum peak energy, E Vy \ , 
and the isotropic equivalent radiated energy, E- lsa , discovered by Amati et al. in 2002 
and confirmed/extended by subsequent osbervations, is one of the most intriguing 
and debated observational evidences in Gamma-Ray Bursts (GRB) astrophysics. In 
this paper I provide an update and a re-analysis of the E p ^ - E- lso correlation basing 
on an updated sample consisting of 41 long GRBs/XRFs with hrm estimates of z 
and observed peak energy, £J p>0 bs > 12 GRBs with uncertain valeus of z and/or E p ^ ODS 
, 2 short GRBs with firm estimates of z and -E Pj0 bs and the peculiar sub-energetic 
events GRB980425/SN1998bw and GRB031203/SN20031w. In addition to standard 
correlation analysis and power-law fitting, the data analysis here reported includes 
a modelization which accounts for sample variance. All 53 classical long GRBs and 
XRFs, including 11 Swift events with published spectral parameters and flucnccs, have 
E Pt i and Ei so values, or upper/lower limits, consistent with the correlation, which 
shows a chance probability as low as ~7xl0 -15 , a slope of ~0.57 (~0.5 when fitting 
by accounting for sample variance) and an extra-Poissonian logarithmic dispersion of 
~0.15, it extends over ^5 orders of magnitude in E- lso and ^3 orders of magnitude in 
E Pt i and holds from the closer to the higher z GRBs. Sub-energetic GRBs (980425 
and possibly 031203) and short GRBs are found to be inconsistent with the E p ^ - 
Ei SO correlation, showing that it can be a powerful tool for discriminating different 
classes of GRBs and understanding their nature and differences. I also discuss the 
main implications of the updated E Py i - Ei so correlation for the models of the physics 
and geometry of GRB emission, its use as a pseudo-redshift estimator and the tests 
of possible selection effects with GRBs of unknown redshift. 
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1 INTRODUCTION 

Since 1997, with the first discoveries of optical counter- 
parts and host galaxies, redshift estimates for Gamma- 
Ray Bursts (GRBs) have become available, allowing the 
study of the intrinsic properties of this challenging astro- 
physical phenomena. Among these, the correlation between 
the photon energy (commonly called peak energy) at which 
the cosmological rest frame v¥ v spectrum peaks, E Pi i , 
and the total isotropic-equivalent radiated energy, _Ei so , is 
one of the most i ntriguing and discu ssed. This correlation 
was discovered bv lAmati et al] J2002T) based on BeppoSAX 
data and subsequently confirmed and extended to X-Ray 
Rich GRBs (XR Rs) and X-Ray Flashes (XRFs) based 
on HETE-2 data jAmatil2003l:lLamb. Donaghv fc Graziaml 
12004 ISakamoto et al ]|2004 l2005alT It can be used to con- 
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strain the parameters of the various scenarios for the physics 
of GRB prompt emission, it is a challenging test for jet and 
GRB/XRF (X-Ray Flashes) unification models and it can 
provide hints on the nature of different sub-classes of GRBs 
(sub-energetic GRBs, short GRBs). Also, the E Pt i - Ei so 
correlation has been used for building up redshift estima- 
tors and is frequently assumed as an imput or as a required 
output for GRB population synthesis models. In this paper, 
I provide an update (up to December 2005) and a re-analysis 
of the E Pt i - Ei so correlation based on a sample of a total of 
56 events which includes Swift GRBs with known redhsift 
and published spectral parameters and two very recent short 
GRBs with firm estimates of redshift and E v ,i ■ The analy- 
sis here reported includes also fitting the data with a model 
which accounts for sample variance, given that this correla- 
tion is highly significant but also shows a dispersion which 
cannot be accounted for only by statistical fluctuations and 
is an important source of information. I also discuss the vari- 
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ous explanations and implications of the existence and prop- 
erties of the -B p ,i - -Eiso correlation, its possible use for the 
estimate of pseudo-redshifts, and the tests based on GRBs 
with unknown redshift aimed to the evaluation of the im- 
pact of possible selection effects. For this last purpouse, I 
also make use of published spectral parameters and fluences 
of a sample of 46 HETE-2 GRBs. 

In order to introduce some basic information for the dis- 
cussion reported in Section 5 and to allow a comparison 
between the results here reported and those reported in pre- 
vious works, the description and properties of the updated 
sample (Section 3) the description and results of the data 
analysis results (Section 4) the discussion of the main im- 
plications and explanations of the E p ,i - E[ so correlation 
(Section 5) and the discussion of pseudo-redshift estimates 
and tests based on GRBs with unknown redshift (Section 6) 
are preceeded (Section 2) by a brief review of spectral and 
energetics properties of GRBs and of previous studies of the 
E p i - -Eiso correlation. 



more than 70 redshift estimates have now become available. 
As a consequence, for these events it is possible to compute 
the intrinsic peak energy E Pi i — E p , b s X (1 + z) and the the 
radiated energy in a given cosmological rest- frame energy 
band fo llowing , e.g., t he m ethods described inlAmati et al 
i 20021) . lAmatil (120031) and iGhirlanda. Ghisellini fc Lazzati 
1 2004). In the simplest assumption of isotropic emission, 
the radiated energy, Ei so , ranges from ~10 50 erg to 
~10 54 erg for most GR Bs and extends down to ^10 49 erg 
when including XRF s llLamb. Donaghv fc Grazianil 120041 : 
ISakamoto et alj|2004l) . When assuming that the GRB emis- 
sion is jet-like, based on achromatic breaks observed in the 
afterglow decay curves of several GRBs, the distribution 
of the collimation corrected radiated energy, E ~ , was ini- 
tially found to be clustered around ~ 10 51 erg (Frai l et alJ 
| 200li iBloom. Frail fc Kulkarnil l2003b : however, recently 
IGhirlanda. Ghisellini fc Lazzatil 11200411 showed that, when 
considering a larger sample of GRBs with known redshift, 



2 THE Ep,i - E iso CORRELATION 

2.1 GRBs peak energy and radiated energy 

The prompt emission spectra of GRBs are non thermal 
and in general c an be modeled with the Band function 
jBand et alJll993t) . a smoothly broken power-law whose pa- 
rameters are the low energy spectral index, a, the high en- 
ergy spectral index, /3, the break energy, Eo, and the over- 
all normalization. In this model, if ft < —2 then the v¥ v 
spectrum peak energy is given by E p , bs = Eo x (2 + 
a). The spectral shape of most GRBs can be satisfactorily 
reproduced by Synchrotron Shock Models (SSM) (e.g. Ta- 
vani 1996 ): the kinetic energy of an ultra-relativistic fire- 
ball (a plasma made of pairs, photons and a small quan- 
tity of baryons) is dissipated into electromagnetic radia- 
tion by means of synchrotron emission originated in in- 
ternal shocks between colliding shells and/or t he exter- 
nal sh ock o f the f i reball with the ISM, see, e.g., iMeszarol 
(2002) and iPiranl d2005l) for recent reviews. Nevertheless, 
the time resolved analysis of BATSE and BeppoSAX GRBs 
showed that, at least during the initial phase of the emis- 
sion, other mechanisms, like Compton up-scattering of UV 
photons surrounding the GRB source by the ultra rel- 
ativistic electrons of the fireball or thermal emission by 
the photosph ere of the fireball, may play an importan t 
role, see, e.g..lPreece etaH ll2000l); iFrontera et al.l feOOOblk 
IGhirlanda. Celotti fc Ghisellinil ll2003t) . The latter emission 
mechanism could also be responsible for the smooth cur- 
vature characterizing GRB average spectra and, in par- 
ticular, may determin e the value of -E p ,obs llRvdd 120051 : 
iRees fc Meszafosl I2005T) . A relevant outcome of the analy- 
sis of BATSE events was the evidence of a substantial clus- 
tering of E Pj0 bs values around 200 keV, but in the recent 
years, the discovery and study of X-ray rich events and X- 
Ray Flashes (XRFs) by BeppoSAX and HETE-2, showed 
that the distribution of E Pl0 bs is much less clustered than 
inferred basing on BATSE data and, in particular, that it is 
characteriz ed by a low energy tail extending down at least 
to ~1 keV dKippen et alJl200lt ISakamoto et al.ll2005bl) . 

Since the BeppoSAX breakthrough discoveries in 1997, 



2.2 Discovery, confirmation and extension of the 
E Pi i — Ei so correlation 

In 2002, I Amati et"aH ||2002|) presented the results of the 
analysis of the average 2-700 keV spectra of 12 BeppoSAX 
GRBs with known redshift (9 firm measurements and 3 
possible values). The more relevant outcome of this work 
was the evidence of a strong correlation between E P; i and 
-Eiso ■ The linear correlation coefficient between log(E Pi i ) 
and log(Eiso ) was found to be 0.949 for the 9 GRBs with 
firm redshift estimates, corresponding to a chance probabil- 
ity of ~0.005%. The slope of the power-law best describ- 
ing the trend of E P; i as a func t ion o f Ei ao was ~0.5. This 
work was extended by lAmatil J2003I) by including in the 
sample 10 more events with known redshift for which new 
spectral data (BeppoSAX events) or published best fit spec- 
tral parameters (BATSE and HETE-2 events) were avail- 
able. The E Pj i - Eiso correlation was confirmed and its sig- 
nificance increased, giving a cor r elation coefficient similar 
to that derived by lAmati et alJ (I2002T) but with a much 
higher number of events. Bas i ng on HETE-2 measurements , 
lLamb, Dona ghy fc Graziani| ( 2004) , ISakamoto etaH (120041) 
and ISak^moto^t^n ~ (l2005a) not only confirmed the E Pi i 

- -Eiso correlation but remarkably extended it to XRFs, 
showing that it holds over three orders of magnitude in 
E p ,i and five orders of magnitude in Ei so • The addition 
of new data, as more redshift estimates became available, 
confirmed the correlation and increased its sign i ficanc e, as 
found e.g. by IGhirlanda. Ghisellini fc Lazzatil <2004ft (29 
events, chance probability of 7.6xl0~ 7 ). Finally, the rel- 
evance of the E Pi i - Eiso correlation for the GRB field 
stimulated several similar studies, which led to the dis- 
covery of correlations of E Pj j with other GRB intensity 
indicator s like the average isotropic equivalent luminos- 
ity Liso llLamb. Donaghv fc Grazianil 120041 l2005h and the 
peak isotropic equivalent lum i nosity L D > n lYonetoku et all 
12004 IGhirlanda et alJl2005al) . lLiang. Dai fc Wul (I2004D also 
showed that, at least for a good fraction of events, the E Pi i 

- Liso correlation holds also within GRBs. All these cor- 
relations show the same slope and dipersion of the E Pl i - 
Eiso correlation, and reflect t he tight correlation existing ^ 
between Ei so , Li so and L P! - lso dLamb. Donaghv fc Grazianil 
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sub-energetic event GRB980425 and the only two short GRBs with firm estimates of z and -E Pi0 b s > GRB050709 and 
GRB051221. The uncertainties are at lcr significance. The "Type" column indicates wether the event is a normal long 
GRB (LONG), an X-Ray Flash (XRF) or is sub— energetic (SUB— EN). The "Instruments" column reports the name of 
the experiment(s), or of the satellite(s), that provided the estimates of spectral parameters and fluence (BAT = BATSE, 
SAX = BeppoSAX , HET = HETE-2, KON = Konus, SWI = Swift ). The last two columns report the references for 
the spectral parameters and the references for the values and uncertainties of £?; so , respectively. The "— " sign indicates 
that Ei so was computed specifically for this work (see text). GRBs detected by Swift are marked with an asterisk. 
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llJlanov et all i2005h. (6)|Price et alJ 1 12002 ) . (7) iBarraud et alJ ll20o1 ), (8)ISakamoto et alJ<2005rJ). (9^ ISakamot^ et alJ 
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(26) IFriedman fc Blooml <2005l) . (271 iGhirlanda et al. I il2005hh. (281 iFenimore et al] <2004l) . (29) Official HETE GRB 
page Ihttp:/ /space. mit.edu/HETE/Bursts/ I. (30l lDotv et al] 120051) ( b ' Uncertainties enlarged in order to account for 
significantly different values measured by different instruments (see text). ( c ) Value derived from time resolved spectral 
analysis by weighting each time interval with the ratio between its fluence and the total GRB fluence. A 20% error is 

rnntaprvFit.i vpl v flscnmpH 
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Table 2. E p ^ and E- lso values for GRBs and XRFs with uncertain estimates, or upper / lower limits, of z or E p<oos ■ The 
uncertainties are at la significance, whereas the upper/lower limits are at 90% c.l. The "Type" column indicates wether 
the event is a normal long GRB (LONG), an X-Ray Flash (XRF) or is sub-energetic (SUB-EN). The "Instruments" 
column reports the name of the experiment(s), or of the satellite(s), that provided the estimates or upper/lower limit to 
the spectral parameters and of the fluence (BAT = BATSE, SAX = BeppoSAX , HET = HETE-2, KON = Konus, SWI 
= Swift ). The last two columns report the references for the spectral parameters and the references for the values and 
uncertainties (or upper / lower limits) of £?; so , respectively. The "— " sign indicates that B; so was computed specifically 
for this work (sec text). GRBs detected by Swift are marked with an asterisk. 
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significantly different values measured by different instruments (see text). ' c ' Upper limit inferred from time resolved 
spectral analysis. 



I2004l2005l : IGhirlanda et al]l2005al) . More intriguing, as will 
be discussed in Section 5, are the E Pi i - E-, and E P] i - Ei so - 
tb (tb is the achromatic af terglow light curve break time) 
correl ation s discovered by IGhirlanda. Ghisellini fc Lazzati 
J2004I1 and lLiang fc Zhand 120051^ basing on a limited sam- 
ple of GRBs with known z, E Pi i and tb- Also, evidence of a 
strong correlation between E Pt i , Ei so and th e "high-si gnal" 
time scale To. 45 has been found very recently llFirmani et alJ 
120061) . 

The outcome of the analysis of the E Pt i - Ei so correla- 
tion performed in previous works are summarized in the first 
five lines of Table 3. As it can be seen, the chance probability 
of this correlation is very low and decreases when increasing 
the number of events in the sample considered. Nevertheless, 
the fits with a power-law are always very poor (as indicated 
by the reported xt values) and both the normalization and 
the index vary significantly depending on the sample con- 
sidered. This is an effect of the extra-Poissonian dispersion 
of the correlation, as will be discussed in detail in the next 
Section. 



3 UPDATED SAMPLE AND DISTRIBUTIONS 

OF £ p ,i AND E iso ■ 

The last three lines of Table 3 report the results of the 
analysis that I performed on the most updated (as of 
December 2005) sample of long GRBs/XRFs with firm 
estimates of both z and E PtOOS . This sample, reported 
in Table 1, consists of a total of 41 events and includes 



events already considered in previous works, new events, 
such as Swift GRBs, and older events for which useful 
spectral information has become available only recently, as 
is the case e.g. for some events detected by Konus/ Wind 
(Ul anov et alJ 120051) . Table 1 includes also the peculiar 
sub-energetic avent GRB980425 and the only two short 
GRBs (GRB050709 and GRB051221) for which firm 
estimates of z and -E p , bs are available. Table 2 includes 
12 classical long GRBs and the other sub-energetic event 
GRB031203, with uncertain estimates of z or E p ^ ODB . It 
also includes XRF030723, for which only an upper limit 



(2002) 




Amati 


ll2003l) IGhirlanda. Ghisellini & Lazzati 


( 2004) Friedman fc Blooml 


(2005), 


Ulanov et alJ 




2005), 



while for the other events they are taken from specific 
references, mostly GCNs. For HETE-2 GRBs with es- 
timate s of -Ep.obs availabl e from both iBarraud et alJ 
ll2003l ) and ISakamoto et alJ d2005aT) . I used the values 



from Sakamoto et al. (2005a), which are based on joint 
WXM/FREGATE analysis. The only exception to this rule 
is GRB020813, for which a depletion of X-ray photons in 
the WXM energy band was observed by m eans of time 
resolved spectral analysis dSato et alJl2005al) . Some GRBs 
included in the samples of some of the works mentioned 
above, were excluded from the sample of GRBs with firm 
estimates of z and E Pt i because of a too large difference 
between the E p ^ ODB measured by different detectors (e.g. 
GRB011121) or a poorly unconstrained 90% c.l. interval of 
E p , obs (e.g. GRB030323) (see references in Table 2). For 
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all GRBs, the last two columns of Tables 1 and 2 report 
the references for the spectral parameters and for the 
values and uncertainties of E p ,i and Ei so , including also 
the methods adopted to derive these quantities based on 
the measured spectra and fluences. It has to be noted that 
in the computation of Ei so some authors assume Ho=65 
km s -1 Mpc -1 (e.g., Amati et al. 2002), while others 
assume Hq=70 km s _1 Mpc -1 (e.g., Ghirlanda, Ghisellini 
& Lazzati 2004); instead, the values of Q. m and Q.a are 
always assumed to be 0.3 and 0.7, respectively. Thus, in 
order to have a homogeneous data set, I corrected the Ei so 
values computed by assuming H =70 km s" 1 Mpc -1 by 
accounting for the different luminosity distance obtained 
with i/ =65 km s" 1 Mpc~\ 

The -Eiso values from lUlanov et all (I2005T) where com- 
puted by integrating the observed fluence between 10/(l+z) 
and 10000/(l+z) keV (Ulanov, private communication) and 
thus between 10-10000 keV in the GRB cosmological rest- 
frame . The difference with respect to integrating from 1 to 
10000 keV clearly depends on the spectral parameters, in 
particular on a, and is not higher than 10% in the worst 
cases (and typically of the ord er of 3-5%). Given t hat the 
values of a axe not reported bv lUlanov et al.| (l2005Tl . I cor- 
rected their Ei so values and la confidence intervals by con- 
servatively assuming that the extension of the integration 
from 10 keV to 1 keV my contribute from to 10% to the 
total radiated energy. 

For some events of the sample, marked with a note in 
Tables 1 and 2, I enlarged the uncertainties on E P: i and/or 
Ei so in order to include the different central values and la 
uncertainties reported in the literature based on measure- 
ments by different GRB detectors. The issue of the sys- 
tematics in the estimates of spectral parameters and flu- 
ence due to detectors limited energy bands and sensitivities 
as a function of photon energy is discussed in Sections 5.3 
and 6.2 . For those few events for which there is no ref- 
erence directly reporting the values of E p ,i and/or Ei so , 
these have been calculated based on published spectral pa- 
rameters and fluence by following the methods mentioned 
in p revious Section and detailed e.g. in jAmati et alJ (2002) 
and iGhirlanda. Ghisellini fc Lazzatil <l2004t) . In particular, 
the unpublished values of -Eiso have been computed in the 
1-10000 keV cosmological rest-frame energy band and by 
assuming a standard cosmology with Sl m =0.3, f2A=0.7 and 
Jfo=65 km s _1 Mpc -1 . For these events, the uncertainties 
on E Pi i and Ei so were computed by propagating the frac- 
ti onal errors on Eo and fluence, respectively, as done, e.g., 
in lAmatil (§003) for BATSE and HETE-2 events. The values 
of z have been taken directly from J. Greiner's GRB Table 1 , 
which also includes complete references. For those cases in 
which only a range for z is available (Table 2), E Pi i and _Ei so 
have been computed by assuming the central value. In the 
case of XR F030723 (Tab l e 2), t he upper limit to the redshift 
comes from |Fvnbo^^dJj 2004) an d the upper lim its to E Pi i 
and Ei so from lLamb. Donaehv fc Grazianil i2004t) . 

It is important to notice that GRBs included in these 
samples have been detected, and their spectral parameters 
measured, by detectors with different sensitivities and en- 
ergy bands. This, together with the fact that Swift is allow- 

1 http:/ /www. mpc. mpg.dc/~jcg/grbgcii. html 



ing z determinations for more types of GRBs, should reduce 
significantly the possible impact of selection effects . This 
issue will be discussed in Section 6. 

In Figure 1 I show the logarithmic distributions of .Eiso (left 
panel) and E Pt i (right panel). These refer to all the GRBs in- 
cluded in Table 1 plus GRB031203 (Table 2). As can be seen, 
in both cases the bulk of the distribution can be fitted by a 
Gaussian, but a low energy tail is evident. In the log(Ei so ) 
distribution, the low tail is due to the sub-energetic GRBs 
980425 and 031203, the XRF 020903 and the short GRB 
050709; the fit of this distribution with a Gaussian gives an 
average of ~10 53 erg and a logarithmic dispersion of ~0.9 . 
It is noticeable that the Ei ao distribution spans about 6 or- 
ders of magnitude. The fit of the log(E Pi i ) distribution with 
a Gaussian gives an average of ~350 keV and a logarithmic 
dispersion of ~0.45 . In this case the low-energy tail is given 
by the XRFs 020903 and 050416a, while sub-energetic and 
short GRBs show E p ,i values consistent with the bulk of the 
distribution. It is worthing to note that the E Pi i distribution 
is much broade r than the E D , bs d istribution inferred from 
BATSE GRBs jKaneko et allteoOrJ) . 



4 THE E p>i - .Eso CORRELATION: 
RE-ANALYSIS 

The (E p ,i ,Eiso ) points corresponding to the 41 GRBs/XRFs 
with with firm estimates of z and E Pi i , all included in Ta- 
ble 1, are shown in Figure 2, whereas the (E p ,i ,Ei so ) points 
and upper/lower limits corresponding to the GRBs with un- 
certain 2 and E Pi i (Table 2) are shown in Figure 3, which 
also includes the points corresponding to the peculiar sub- 
energetic GRB980425 and the two short GRBs 050709 and 
051221 . In both Figures, the point corresponding to Swift 
GRBs are shown as filled circles. The first two lines of the 
second part of Table 3 report the results of the analysis per- 
formed on the sample plotted in Figure 2. The correlation 
analysis is based on the estimate of the Spearman's rank cor- 
relation coefficient between E p ,i and _Ei so and the fits with 
a power-law E p .i = K x E™ are performed by account- 
ing for the errors on both E p ,i and Ei so • As can be seen, 
with the updated sample subject of this analysis, the chance 
probability of correlation between the logarithms of E p ,i and 
Eiso is as low as ~3.5xl0 -14 , when considering only the 41 
classical GRBs, and ~10 -14 when including also the XRFs 
020903 and 050416a. Thus, increasing the sample by adding 
new data and making it more complete (as mentioned above) 
not only confirms the E Pi i - Ei so correlation but also reduces 
its chance probability to a negligible value. The index of the 
power-law, ~0.57, is found to lie in the range 0.4-0.6, con- 
sistently with the findings of previous analysis, and does not 
change significantly by including or not in the sample XRFs 
020903 and 050416a (see Table 3). The value of the nor- 
malization is found to be somewhat lower with respect to 
previous analysis, except for the recent analysis performed 
by Nava et al. (2006) on a subsample of 18 GRBs which 
includes more recent events with respect to previous works. 
The power-law best fitting the data of the 41 GRBs with 
firm estimates of z and E p ,i is shown as a dashed line in 
Figures 2 and 3. 

Despite the correlation is very highly significant, the x 2 
values obtained by fitting the data with a power-law are 
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Figure 1. Distributions of log(i?i so ) (left) and log(_E Pi i ) (right) for the 41 GRBs with firm redshift and -B Pi0 bs > the short GRBs 050709 
and 051221 and the sub-energetic GRBs 980425 and 031203. For both distributions, the best fit Gaussian is superimposed to the data. 



Table 3. Summary of the results of the analysis of the E p \ - £j so correlation as reported in previous works (top) and 
performed in this work (bottom). The coefficient p is the Spearman's rank correlation coefficient between E Pt i and E, so . 
N is the number of events considered, m, K and x 2 refer to fits of the E p j - Ej so correlation with a power-law £ Pi ; = 
KxE™ (E P: i is in keV and i?i so in units of 10 52 erg). The uncertainties reported in the first part of the Table (values 
taken from the literature) are at 1<t confidence level, while those reported in the second part (results obtained in this 
work) are at 90% significance. When not available in the literature, the values of K and x 2 reported in the top part of 
the Table have been computed specifically for this work. 



Reference 


N 


P 


Chance Prob. 


m 


K 




Amati et al. (2002) 
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0.92 


5.0X10" 4 


0.52±0.06 


105±11 


3.9 


Amati (2003) 


20 


0.92 


l.lxlO" 8 


0.35±0.06 


118±9 


6.1 


Ghirlanda, Ghisellini & Lazzati (2004) 


27 


0.80 


7.6xl0" 7 


0.40±0.05 
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Friedman & Bloom (2005) 
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0.88 


4.9xl0" 10 
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Nava et al. (2006) 


IS 


0.82 
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5.2 


This work (only GRBs) 

This work (including XRFs 020903 and 050416) 
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poor, as found in previous analysis (Table 3). This means 
that the scatter of the data around the best fit model cannot 
be due only to statistical fluctuations, unless the systemat- 
ics in the estimates of E p ,i and Ei so are strongly under- 
estimated, as will be discussed in next Section. Another ef- 
fect of the dispersion characterizing the correlation is that, 
as mentioned above, the slope and normalization of the 
power-law are found to change significantly depending on 
the sub-sample considered. This extra-Poissonian disper- 
sion of the E Pt i - Ei so correlation, which potentially contains 
precious information (as will be discussed in Section 5) and 
has to be taken into account when testing it (as will be dis- 
cussed in Section 6), can be quantified by introducing in the 
modelization a further parameter commonly called "sample 
variance" or "slop". The issue of fitting data with a power- 
law by accounting simultaneously for X an d Y errors and 
for sa mple variance has been faced e.g. b y iReichart et alJ 
feOOll) and iReicharTfc Nvsewanderl J2006T) when analyzing 
the peak luminosity-variability correlation in GRBs. 

The methods used in these works are based on a likeli- 
hood function derived with a bayesian approach to the prob- 



lem; h owever, recentlv |D ' Agostinil (I2005T) and lGuidorzi et ail 
(2006) showed that the correct likelihood function i s 
slight ly different from that use d by IReichart et al] i200ll) ; 
Reichart fc N Y sewan der| J2006T). I applied t he method by 
D'Aeostinil | |200E|) and lGuidorzi et all J200fj) to the sample 
of 41 GRBs considered above; with this modelization the pa- 
rameters are the index and normalization of the power-law 
(m and K) and the logarithmic dispersion of E Pi i (ai og E p ()■ 
The result of this analysis is reported in the last line of Table 
3. The values of the index and normalization of the best-fit 
power-law, ~0.5 and ~100, respectively, lie in the ranges of 
values found with different samples by adopting the simple 
power-law fit and are coincident with those usually assumed 
in the literature when comparing new data with the E Pt i - 
Ei ao correlation or using it as an input or required output 
for GRBs/XRFs synthes is models (see next Section), basing 
on the early results from I Amati et ail J2002) ( Table 1, first 
line). The best fit power-law obtained with this method is 
plotted in Figures 2 and 3 as a continuous line. The value 
of the sample variance resulting from the fit is ai og E p ( = 
°- 15 ±o.04 ( 90% °- L ); in Figures 1 and 2 I show the region, 
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delimitated by two dotted lines, corresponding to deviations 
of E Pt i from the best fit power-law of ~2.5ai og E p f , assuming 
o~iogE p s =0.15 (the central value of the 90% confidence inter- 
val) . For a comparison, as shown in Figure 4, the dispersion 
of the log(_E Pi i ) central values around the best fit power- 
law obtained without accounting for sample variance (line 
7 of Table 3) can be fitted with a Gaussian with dispersion 
~0.21, consistently with previous analysis based on smaller 
samples. A similar value {(Ji og E 4 ~0.2, is obatined when 
computing the scatter of the central data points around the 
law _E p .i = 95 x Ef S Q 9 . This value is of course higher than 
sample variance, because it includes statistical fluctuations; 
if one considers this as the overall dispersion of the correla- 
tion, then the two dotted lines in Figures 2 and 3 delimitate 
the ±2<r region. An estimate of the sample variance char- 
acterizing the E u .j - Ej ao correlatio n was also performed by 
lLamb. Donaehv fc Grazianil ({2005) , who, based on a smaller 
sample and following a different method, derived a value of 
-0.13 . 

Finally, from Figure 3, it can be seen that the uncertain 
values and upper/lower limits of E Pi i and Ei ao of classical 
long GRBs and XRFs (reported in Table 2) are consistent 
with the E P j - Ei BO correlation, including Swift GRBs. Fig- 
ure 3 also shows clearly that the E Pi i - Ei so plane can be 
very useful in discriminating different classes of GRBs. In- 
deed, both sub-energetic GRBs (GRB980425 and, possibly, 
GRB031203) and short GRBs (050709 and 051221) are clear 
outliers of the correlations, showing i?i so values too low with 
respect to their peak energies, which range within those of 
normal GRBs (Figure 1, right panel). 

During the reviewing process of this article, firm es- 
timates of both z and -E P!0 b s have become available for 4 
more Swift events (GRB060115, GRB060124, GRB060206, 
GRB060418), based on spectral measurements from Swift, 
HETE-2, Konus. The E Pt i and _Bi so of these events are all 
fully consistent with the E p ,i - Ei so correlation. 



5 MAIN IMPLICATIONS AND DISCUSSION 

The analysis presented in previous Sections, based on an 
updated sample containing about twice events with re- 
spect to previous works and including the recent Swift 
GRBs, confirms and strenghtens the E Pi i - E[ so correla- 
tion for long GRBs/XRFs and gives its best characteri- 
zation up to now in terms of index and normalization of 
the best-fit power-law and of its dispersion. Remarkably, 
the correlation now extends over ~5 orders of magnitude 
in Ei so , ~3 orders of magnitude in E p ^ and over a red- 
shift range ~0.15<z<4.5 (but also the highest z event, 
GRB050904 at z—6.29, has E p i and -Ei so values consis- 



tent with it). Since its discovery in 2002 dAmati et alJ 
2002) and in particular its confirmation and extension 
to XRFs |Amatill2003l : lLamb. Donaehv fc Grazianil 12004 
ISakamoto et alJl2004F 7 the origin of the E Pi \ - E\ so corre- 
lation and its implications for GRB models have been in- 
vestigated by several works. The impact of this robust ob- 
servational evidence on prompt emission models concerns 
mainly the physics, the geometry (i.e. shape and proper- 
ties of jets), viewing angle effects and GRB/XRF unifica- 
tion. Indeed, the existence of the _E Pj j - _Bi so correlation and 
its properties are also often used as an ingredient or a test 



output for GRB synthesi s models, as in the case of, e.g., 
the GRB/XRF model bvlBarraud et all (120051) . the m ulti- 
subjets model bv lTomaT Yamazaki & Nakamura ( 200 51), th e 
uniform jet model by lLamb. Donaehv fc Grazianil (2005), 
the st udy of the imp act of off-jet relativistic kinemat- 
ics by iDonaehvl d2006l) . the dissipative photosphere mod- 
els bvj^e^^^^eizafogj20o3) i _the supercritical pile model 
bv lMasticmadi^^^^azanasI (1200(11) . Another important out- 
come of the analysis presented in this paper is the clear 
evidence that, in addition to the peculiar sub-energetic 
GRB980425 (and possibly GRB031203), short GRBs do not 
follow the E Pi i - Ei so correlation, as suggested by the differ- 
ent location between long and short BATSE GRBs in the 
hardness-intensity plane. Below I summarize these topics 
and discuss also the possible origin of the extra-Poissonian 
dispersion of the correlation. 



5.1 Physics of prompt emission 

The physics of the prompt emission of GRBs is still far to 
be settled and a variety of scenarios, within the standard 
fireball picture, have been proposed, based on different emis- 
sion mechanisms (e.g. SSM internal shocks, Inverse Comp- 
ton dominated internal shocks, SSM external shocks, pho- 
tospheric emission dominated models) and different kinds 
of fireball (e.g. kine tic energy dominat e d or Poynting flux 
dominated), see e.e. lZhane fc Mcszaros d2002l) for a review. 
In general, both E Pt i and -Eiso are finked to the fireball bulk 
Lorentz factor, T, in a way that varies in each scenario, and 
the existence and properties of the E Pt i - Ei so correlation al- 
low to constrain the range of value s of the parame ters, see, 
e.e.. lZhane fc Meszarosl d2002Fl and ISchaefej j|2003h For in- 
stance , as shown, e.e.. bv |^angfc Mesz^ro^ (12002! ) . fRvdd 
d2005l) and iRees fc Meszarosl j200^T " for a power-law elec- 
tron distribution generated in an internal shock within a 
fireball with bulk Lorentz factor V, it is possible to derive the 
relation E p ,i oc r~ 2 L 1,/2 f~ 1 , where L is the GRB luminosity 
and t v is the typical variability time scale. Clearly, in order 



to produce the observed i? Pi i - E\ so correlation the above 
formula would require that V and t v are approximately the 
same for all GRBs, an assumption which is difficult to jus- 
tify. Things get even more complicated if one takes into ac- 
count that the models generally assume L oc T 13 , with the 
va lue of varying in each scenario and is typically ~2- 
3 dZhane fc Meszarosl [20021: ISchaefeJ l2003t iRamirez-RuTzl 
l2^05[lRvd^l2005l) . More specific examples of the constraints 
put by the E P:1 - Ei ao correlation on the parameters of 
SSM and IC based emission models, b oth in internal and 
extern a l shocks, can b e found, e.g., in IZhane fc Meszarosl 
(2002); Schacfei] d2003l) . An interesting possibility, which 
is currently the subject of many theoretical works, is that 
a substantial contribution to prompt radiation of GRBs 
comes from direct or Comptoniz ed thermal emission from 
the photosphere of the fireba ll ||ZIuuig_fc_Mes^aroj |2002j; 
iRvdel 120051: IRees fc Meszarosl 120051: iRamirez-RuizI booslf . 
This could explain the ve ry hard spectra observed at the be- 
ginnin g of several evente jjPreece^^lJ|20poyR ontera et aid 
boOObl: iGhirlanda. Celotti fc Ghisellind 120031) . inconsistent 
with SSM models, and the smooth curvature characteriz- 
ing GRBs average spectra. In this scenarios, E Pt i is mainly 
determined by the peak temperature T p k of black-body dis- 
tributed photons and thus naturally linked to the mini- 
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Figure 2. E Pil and Ei so values for 41 GRBs/XRFs with firm redshift and -E Pi0 b s estimates. Filled circles correspond to Swift GRBs. 
The continuous line is the best fit power— law E p \ = 95 X _EP' 49 obtained by accounting for sample variance; the dotted lines delimitate 
the region corresponding to a vertical logarithmic deviation of 0.4. The dashed line is the best fit power-law E p j = 77 X EPj? 7 obtained 
by fitting the data without accounting for sample variance. See text for details. 




Figure 3. Same as Figure 2 for the 12 GRBs with uncertain estimates of z and/or -E Pj0 bs , for the peculiar sub— energetic event 
GRB980425 / SN1998bw and for the two short GRBs 050709 and 051221. 



nosity or radiated energy. For instance, for Comptonized 
emission from the photosphere one can derive the relations 
E p ,i oc FT pk oc r 2 L" 1/4 or E p ,i oc TT pk oc r<T 1/2 L 1/4 
(where ro is a particular distance between the central en- 
gine and th e emitting region] . dep ending on the assump- 
tions made jRees fc Meszafosll2005l) . Also in this case, off 
course, the L oc T 13 relation plays a decisive role. As shown 



bv lRees fc Meszafosl 1)20051) . in this scenario the correct E p s 
- Ei BO relation can be obtained for some specific physical 
conditions just below the photosphere. 

Finally, also, the fact that the E Pi i distribution is 
broader than inferred previously basing mainly on the ob- 
served -E Pj0 bs values of bright BATSE GRBs, as shown in 
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Figure 4. Logarithmic dispersion of the E Pi i values of the 41 
GRBs with firm redshift and E p j estimates around the power- 
law best fitting the E Ptl - Ei BO correlation without accounting for 
sample variance. 1 also show the best fit Gaussian, which has a 
dispersion of cr log ( Ep i) ~ 0.21 . 



Figure 1, can put important con straints on models for th e 
physics of GRB prompt emission JZhang & Mcszaros 2002). 



5.2 Jets, viewing angle effects and GRB/XRF 
unification models 

The validity of the E Ptl - E- lso correlation from the most 
energetic GRBs to XRFs (see Figure 2) confirms that these 
two phenomena have the same origin and is a very chal- 
lenging observable for GRB jet models. Indeed, these mod- 
els have to explain not only how Eiso and E Pt i are linked 
to the jet opening angle, 9j e t, and/or to the viewing an- 
gle with respect to the jet axis, 8 V , but also how Ei so 
can span over several orders of magnitu des. In the most 
simpl e scenario, the uniform jet model jFrail et alJl200lt 
lLamb. Donaghv fc Grazianil l2005t) . jet opening angles are 
variable and the observer measures the same value of Ei so 
independently of 9 V . In the other popular scenario, the 
universal structured jet model (e.g. Rossi, Lazzati & Rees 
), E\ so depends on 6 V . As discussed in Section 2, in the 
hypothesis that achromatic breaks found in the afterglow 
light curves of some GRBs with known redsh ift are due 
to co l limated emission, it was origin ally found dFrail et alJ 
1200 it iBereer. Kulkarni fc Frail I2003T) that the collimation 
corrected radiated energy, Ej , is of the same order (~10 51 
erg) for most GRBs and that Eiso oc 0j^ t , assuming a uni- 
form jet. In the case of structured jet models, which assume 
that 6 j e t is similar for all GRBs (hence this scenario is also 
called universal jet model) the same observations imply that 
i?iso oc 8~ 2 . Thus, always under the assumption of a nearly 
constant Z? 7 , the found E Pt i - Ei so correlation implies E Pi i oc 
and -E Pj j oc 9~* for the uniform and structured jet mod- 



"jet "P,l ^ " V .,^ upu^i, 

els, respectively. Lamb, Donaghv & Graziani (2005) argue 
that the structured universal jet model, in order to explain 
the validity of the E Pt i - Ei ao correlation from XRFs to en- 
ergetic GRBs, predicts a number of detected XRFs several 
orders of magnitude higher than the observed one (~l/3 
than that of GRBs). In their view, the uniform jet model 



can overcome these problems by assuming a distribution of 
jet opening angles N(6j e t) oc 9J^ t . This implies that the great 
majority of GRBs have opening angles smaller than ~1° and 
that the true rate of GRBs is several orders of magnitude 
higher than obs erved and comparab le to that of SN Ic. On 
the other hand, IZhang et alJ <l2004l) show that the require- 
ment that most GRBs have jet opening angles less than 1 
degree, needed in the uniform jet scenario in order to ex- 
plain the E Pi i - Ei ao correlation, as discussed above, implies 
values of the fireball kinetic energy and/or of the interstellar 
medium density much higher than those inferred from the af- 
terglow decay light curves. Together w ith other auth o rs, e.g. . 
iLlovd-Ronning. Dai fc Zhand (120041) : iDai fc Zhand" <2003h 
they propose a modification of the universal structured jet 
model, the quasi-universal Gaussian structured jet. In this 
model, the measured Ei so undergoes a mild variation for val- 
ues of V inside a typical angle, which has a quasi-universal 
value for all GRBs/XRFs, whereas it decreases very rapidly 
(e.g. exponentially) for values outside the typical angle. In 
this way, the universal structured jet scenario can reproduce 
the E Pt i - Ei S o correlation and predict the observed ratio be- 
tween the number of XRFs and that of GRBs. Recently, a 
Fisher-shape has been propose d, for both the variable angle 
and universal angle scen arios iDonaghv. Graziani fc Lambl 
120051: IDai fc Zhand" [2005 ) , as a very promising alternative, 
in particular for the explanation of the validity of the E P: i — 
Ei S o correlation from the brightest GRBs to XRFs. Other jet 
models proposed very recently that can reproduce the E Pt i 
- Ejso correlation include the ring-shaped jet model, see 
lEichler fc Levin son| I^O O"^. and the multi-component (sub- 
lets') model. see lTmiia" Yamazaki fc NakamuraN2005l) . 

Of particular interest are the off-axis scenarios, in which 
the jet is typically assumed to be uniform but, due to rel- 
ativistic beaming and Doppler effects, for 9 V > 9 jet the 
measured emissivity does not sharply go to zero and the 
event is detected by the obse rver with Ej B o and E p ,j drop 



,3ing rapidly as 9 V increases \ Yam azaki. Ioka & Nakamura 
| 2003[lGranot et al.l2002l : lEichler fc Levinsonl2004l : lDonaghv 

2006). In these models, XRFs are those events seen very 
off-axis and the XRFs rate with respect to GRBs and 
the E p ,i - Ejso correlation can be correct l y pre dicted. As 
shown e.g. in I Yamazaki. Ioka fc Nakamural (2004) for a sim- 
ple model of GRB jet, if the Doppler shift factor is S — 
[r(l — f3cos(9 v — 9 jet)]' 1 (where j3 is the velocity of the out- 
flow in units of speed of light), E Pt i and Ei so scale, with re- 
spect to their values observable at the edge of the jet, as E p s 
oc 8 and Ei so oc S 1 '" , where a is the spectral index of the 
prompt emission photon spectrum in the hard X-ray energy 
band. By combining these relations one can obtain the E Pt i 

- Eiso correlation with index 0.5 for classical GRBs (q~— 1) 
and 0.3 for XRFs (a~— 2). A detailed study of off-jet rel- 
ativistic kinema tics effects has been recently performed by 
iDonaghvl ll200fj) for a set uniform (i.e. top hat shaped - vari- 
able opening angle) jet models, finding that these scenarios 
predict a significant population of bursts away from the E Pt i 

- Eiso correlation, unless T > 300 for all bursts or there is 
a strong anti-correlation between F and the jet solid angle. 
Finally, the off-axis effects for very weak and soft events 
can be applied in a similar way as described above in the 
context of the cannon ball (CB) model for G RBs in order 
to re produce the E Pl i - Ei ao correlation jDar fc De Ruiulal 
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5.3 The dispersion of the correlation 

In addition to its existence and slope, also the extra- 
Poissonian dispersion of the E Pi i - Ei so correlation is an 
important source of information. As discussed in Section 3 
and shown in Table 3, while the correlation is very highly 
significant, the scatter of the data around the best fit power- 
law exceed that expected by statistical fluctuations alone 
and produces high values of xt ■ By fitting with a Gaussian 
the dispersion of the central values of log(E p .i ) around the 
best fit model, I obtain <Ji og E p ,- ~0-2, while by fit ting the 
whole data with the method by iD'Agostinil i2005l) , which 
includes sample variance directly in the model, I obtain 
ciogE i = 0.151q 4 . A similar scatter, even if computed 
with only the first of the two methods reported above, is 
found for the E D ,j - Lj BO and E D ,j - L poa k,i so correlations, 
see, e.g.. iGhirlanda et alJ (l2005al) . Intriguingly, the E Pl j - 
E 7 correlation shows instead a lower dispersion, of the or- 
der of ffiooE. .-~ 0-l iGhirlanda. Ghisellini fc Lazzatil 120041: 
iNava et alJl2006D . even if this correlation is based on a still 
low number of events, it requires an estimate of tt in addition 
to i?p,obs and z, it depends on jet model and circum-burst 
environment properties (density, distrib ution) and there are 
possib le outliers, as discussed e.g. by IFriedman fc BloomI 
(2005) and may be indicated by the lack of break in the X- 
ray afterglow light curve of some Swift GRBs with known 
redshift. A low 3-D dispersion characterizes also the E Pi i - 
Eiso - tb correlation, which is a kind of model-independent 
version of the E D ,j - E 1 correlation llLiang fc Zhand 120051: 
INava et alll200d) . The existence of both the E p ,j - Ej ao and 
E Pt i - E 7 correlations is due to the fact that the collima- 
tion angles of GRBs are distributed over a relatively narrow 
range of values; the lower dispersion of the E p , i — E*y corre- 
lation indicates that at least part of the scatter of the E Pi i 
- -Eiao correlation is due to the dispersion of jet opening an- 
gles. And indeed, the comparison of the properties of the two 
correlations has been used, in addition to the study of the 
relation between jet opening angle and radiated energy, to 
inf er the distribution of jet opening angle s, as done, e.g., 
bv lGhir landa 1 Ghisellin i fc Firmanil ll2005l) : iBosniak et al.l 
( 2006a) : iDonaghvl 120061) . Very recently, it has been found 
evidence that a relevant contribution to the dispersion of 
the correlation is due to temporal properties of the prompt 
emiss ion, like the half-width of the auto-correlation func- 
tion iBorgqnovo fc Bj ornsson||2006l) and the "high-signal" 
time scale jFirmani et al.ll2006l) . Other contribution to the 
scatter of the E D , j - Ej ao correlation may com e from viewing 
angle effects, e.g. lLevinson fc Eichlerl i2005l) . the dispersion 
of the parameters of the fireball and/or of the time scales 
(as discussed in Section 5.1 concerning synchrotron emis- 
sion in internal shocks), the inhomogeneous s tructure of the 
jet, e.g. lToma. Yamazaki fc Nakamural <l2005|) ). the possible 
presence of significant amount of material in the circum- 
burst region, that would affect the estimates of both Ei ao 
and E Pt i with a global qualitative effec t of steepening th e 
correlation and increasing its dispersion jLongo et al . 2005). 
In general, several GRB population synthesis models, like 
those mentioned at the beginning of this Section, predict a 
scatter of the E p ,i - Ei so correlation which depend on the 
parameters values. 

When investigating the above physical interpretations 
and implications, it is important to take into account that 



at least part of the dispersion of the E Pi i - E[ so correla- 
tion could arise from instrumental and other systematic 
effects in the estimates of E D ,j and Ej so and their uncer- 
taintie s . As discussed e.g. by jLlovd. Petrosian fc Mallozzi 
( 2000); iLlovd fc Petrosian! (I2002T) . data truncation effects, 
i.e. the systematics introduced by the limited energy band of 
the detector, may affect significantly the estimate of E p , bs 
. Indeed, as discussed in Section 2, typical GRB spectra are 
characterized by a very smooth curvature and cover ~3 or- 
ders of magnitude or more in photon energy. Thus, unless 
the energy band extends from few keV to few MeV, only 
a portion of the spectrum can be detected by a single in- 
strument, which may cause a bias in the estimate of the 
spectral parameters, especially when E PlC ,bs is not far from 
the low or high thresholds. Also, both BeppoSAX (~2-700 
keV) and HETE-2 (2-400 keV), in addition to be capable 
to detect X-ray rich events and XRFs which could not be 
triggered by BATSE (~25-2000 keV), showed that X-ray 
emission below ~30 keV of normal GRBs can last up to 
several tens of seconds more than in the hard X-ray energy 
band. Thus, a GRB detector working at energies higher than 
few tens of keV, like BATSE, may have lost, for a fraction of 
GRBs, a substantial portion of the soft X-ray emission, with 
a consequent overestimate of a and E Pj0 b s • These effects 
can indeed be seen when comparing the X - and hard X- 
rays d u ration and light cur ves of BeppoSAX iFrontera^^S] 
| 2000b l : lAmati et aiTl2002l) and HETE-2 JSakamoto et alJ 
2005b) GRBs and the average spectral paramet ers estimated 
by BeppoSAX /WFC+GRBM and BATSE (lAmati et alJ 
2002; Ijimenez. Band fc Piranl l200lT) for those events re- 
vealed by both satellites. This is true, even if to a mi- 
nor extent, when comparing the best fit spectral mod- 
els obtained with HE TE-2 /FREG ATE (7-400 keV) alone 
fearraud et alj|2003l) with those obtained by jointly fitting 
HETE-2/WXM (2-30 k eV) and HETE-2/FREGATE data 
JSakamoto et alJl2005bh . 

Concerning Ej so , the main source of possible system- 
atics is the extrapolation to the 1-10000 keV cosmological 
rest-frame energy band of the spectral model obtained by 
fitting data in the instrument energy band (see Section 2). 
Indeed, the (known) statistical uncertainties and (unknown) 
biases in the estimates of spectral parameters may affect 
significantly the estimate of Ei so . This is particularly true 
for estimates based on spectra from instruments with high 
energy bound at a few hundreds of keV, like HETE-2 or 
Swift /BAT, which in several cases cannot provide a reli- 
able estimate of the high energy spectral index f3 . In ad- 
dition, the typical choice of computing Ei ao in the 1-10000 
keV rest-frame energy band may not be optimal for very 
soft events with values of E Pi i below a few tens of keV, for 
which this method can likely lead to an underestimate of 
-Ei so • Off course, also the choice of the cosmological param- 
eters for the computation of the luminosity distance, usually 
made by assuming values in the ranges given by the so called 
"concordance cosmology" based on type la SNe and CMB 
measurements, affects the values of Ei so . 

Finally, very recently Swift /XRT found evidence of X- 
ray flares following the en d of the promp t emis sion as de- 
tected by Swift /BAT, e.g. iBurrows etaU (2005). One pos- 
sible interpretation of these phenomena is that they are 
due to con tinued activity o f the engine and /or late inter- 
nal shocks <Fan fc Weill2005UBurrows et alJlioollWu et all 
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l2006tlPerna et aJjl2006tlProga fc Zhandl200d) and thus can 
be considered part of the prompt emission. Given that these 
events are typically soft and that their fluence can be a sig- 
nificant fraction of that of the GRB, their non-detection 
(because of sensitivity) by past and current GRB detectors 
may also have biased the estimates of Ei so (under-estimate) 
and of E Pt i (over-estimate). 

5.4 Short GRBs 

In the past, the E Pt i - Ei so correlation has been studied 
basing on data of long GRBs, given that no redshift infor- 
mation was available for short GRBs. Nevertheless, thanks 
to the measurements performed by HETE-2 and Swift , in 
the last year it has been possible to detect afterglow emis- 
sion and possible optical counterparts and/or host galaxies 
for a few short GRBs. As discussed in Section 3, I have in- 
cluded in the analysis reported in this paper the two short 
GRB s with firm estimates of z and -Ep. obs : GRB050709 
Villasenor etafllgOOSt iHiorth et alj|2005h and GR B051221 
Golenetskii et alJl2005et iBerger fc Soderber J2005I) . As can 
be clearly seen in Figure 3, these events are outliers to 
the _E P; i - Eiso correlation, even when taking into account 
its extra-Poissonian dispersion. In addition, the spectral 
data of other recently l ocalized short GRBs with po ssible 
redsh ifts, GRB0505 09b dGehrels et alJl2005HPedersen et ap 
120051) . GRB05072 4 dKrimm et al]l200St iBereer et alJl2005h 
and GRB050813 l|Sato et alJl2005bl lBergedl200ifl . indicate 
a likely inconsistence with the E Pt j - E- lso correlation, even 
though an estimate of -E p , bs was not possible for these 
events. These evidences confirm the expectations based on 
the fact that short GRBs tend to form a separate class with 
respect to long GRBs in the hardness-intensity plane, i.e. 
they tend to be weaker and harder, and clearly shows the po- 
tentiality of the use of the E Pt i - Ei so plane for distinguishing 
different classes of GRBs and understanding their nature. In 
particular, given that short and long GRBs partially over- 
lap in the hardness-intensity and hardness-duration planes, 
that a fraction of short GRBs show a soft er extended emis- 
sion w hich can last tens of seconds, e.g. iNorris fc Bonnehl 
(2006), the possible relevant impact of spectral/temporal 
trigger selection effects, it is sometimes difficult to estab- 
lish to which class a GRB belongs. Thus the E Pt i - Ei so 
plane may be a powerful tool under this respect, espe- 
cially when combined e.g with the lack of spectral lag and 
spectral evolution which seems to characterize short GRBs 
iNorris fc Bonnellll2006h . 

While the commonly accepted hypothesis for the pro- 
genitors of long GRBs is their association with the core col- 
lapse of massive fastly rotating stars (so called "hypernova" 
or "collapsar" models), based on their duration (from few 
seconds up to hundreds of seconds), their huge _Bi so (Fig- 
ure 1), their typical location inside blue galaxies with high 
star formation rate and the evidence of a metal-rich circum- 
burst environment (as inferred from absorption / emission 
features in prompt and afterglow emission spectra), short 
GRBs are thought to originate from the coalescence of neu- 
tron star - neutron star or neutron star - black hole binaries 
(the "merger" scenarios). It has also been proposed that a 
fraction of them may be gia nt flare from ext r a-galac t ic soft 
gamma repeaters, see, e.g.. iMeszaros! (120021 ); IPiranl (2005) 
for reviews. The very recent observations mentioned above, 



show that, as long GRBs, also short GRBs lie at cosmologi- 
cal distances (0.1<z<l) but they are less energetic (see Ta- 
ble 1 and Figures 1 and 3). Also, their host galaxies show dif- 
ferent morphologies and star forming activity: for instance, 
short GRB050724 (as possibly short GRB050509b) came 
from an elliptical galaxy with low star formation, whereas 
short GRB050709 was associated with an irregular late type 
star forming galaxy. Both events were located towards the 
outskirts of their host galaxies. All these properties match 
the predictions of the merger scenarios (e.g. Belczynski et 
al. 2006 ). Off course, these evidences concern still a very 
low number of short GRBs and thus it is premature to 
draw any definitive conclusion. Anyway, the different nature 
of the progenitors between short and long GRBs can help 
in un derstanding their different behav i or in the E Pl i - Ei BO 
plane. iGhirlanda. Ghisellini fc Celottil d2004) found, basing 
on BATSE data, that the emission properties of short GRBs 
are similar to those of the first ~l-2 s of long GRBs. This 
could indicate that the central engine is the same for the 
two classes, but works for a longer time in long GRBs. This 
would explain the low radiated energy by short GRBs and, 
given that the emission would stop before the typical hard to 
soft evolution observed for long GRBs, also their high E p j 
(with respect to their Ei so ). The merger scenarios naturally 
explain the short life of the central engine, and thus the 
low radiated energy and the quitting of the emission before 
hard to soft spectral evolution. In addition, they a lso predict 
a wea k afterglow emission, as recently observed llFox et alJ 
2005), because of the cleaner and lower density circum-burst 
medium with respect to that predicted by the hypernova sce- 
narios for long GRBs, which would cause a very inefficient 
external shock. In the GRB scenarios where most of the 
prompt emission of long GRBs is due to the external shock 
too, this naturally explains the lack of long lasting and soft- 
ening emission in short GRBs. Again, this would produce 
an high average E p ^ value with respect to the radiate en- 
ergy, and thus the inconsistency with the E Pi i - Ei so corre- 
lation holding for long GRBs. Finally, in order to explain 
the very low or spec tral lag observed in shor t GRBs with 
respect to long GRBs. lNorris fc Bo nncll (2006) consider the 
hypothesis that the typical F of short events is several times 
that of long ones, and thus of the o rder of ^500-1000, as 
predic ted e.g by the merger model of lAlov. Janka fc Mulled 
(2005). With such a high Lorentz factor, internal shocks are 
expected to have a low efficiency, and indeed this is one pos- 
sible explanation for the weakness and softness of XRFs (e.g 
Barrau d et al. 2005 ). If we assume the scen ario proposed 
by e.g. IGhirlanda. Celotti fc Ghisellinil d2003l) . in which the 
spectrally hard emission observed in the first seconds of long 
GRBs is due to (possibly Compton dragged) thermal emis- 
sion from the photosphere and the later emission to syn- 
chrotron processes occurring in internal shocks, the high F 
would thus produce a short-hard emission possibly followed 
by a weak soft component, as observed for several short 
GRBs dVillasenor et alJl2005t INorris fc BonnelJl2006|) . 

5.5 Sub-energetic GRBs and the GRB/SN 
connection 

As can be seen in Figure 3, in addition to the two short 
GRBs also the prototype event for the GRB/SN connec- 
tion, GRB980425/SN1998bw, is characterized by values of 
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E Pt i and Ei so completely inconsistent with the E Pt i - Ei so 
correlation holding for the other events. From an obser- 
vational point of view, this is a direct consequence of the 
fact that the event is characterized by a fluence and a 
measured peak energy in the range of classical long GRBs 
but, based on the commonly accepted association with 
SN1998bw, it lies at a much lower distance (z = 0.0085). 
Figure 3 shows that also another event associated with 
a SN event, GRB031203, is characterized by a value of 
_Ep.i which, combined with its low value of Ei so , makes 
it completely inconsistent with the correlation. Given that 
GRB031203 is the most similar event to GRB980425 un- 
der several points of view (although lying at a larger dis- 
tance, z ~ 0.1), in particular the strong evidence of associ- 
ation with a SN event (SN20031w) and t he low afterglow 
energy inferred from radio observations llSoderberg et alJ 
120041 : ISazonov. Lutovinov fc Sunvaevl l2004f) . this inconsis- 
tency has been invoked as a further evidence of the ex- 
istence of a class of close sub-energetic GRBs. How- 
ever, the lower limit on E p .j bas ed on ISGRI data 
iSazonov. Lutovinov fc Sunvaevl l2004f) and the E Pi i esti- 
mate based on Konus data (Table 2, Ulanov et al. 2005) 
are currently debated based on the dust echo observed 
with XMM-Newton, which could indicate a m uch softer 
prompt emission spectrum dWatson et alJ 1200(f) . The sam- 
ple of GRBs with most evidence of association with a 
SN include also GRB030329 (SN2003dh) and GRB021211 
(SN20021t), which, in converse, are not sub-energetic and 
are characterized by E Pt i and -Eiso values fully consistent 
with the E Pt i - Ej so correlation. The fact the two closest 
and sub-energetic among those GRBs most clearly associ- 
ated with a SN are outliers to the E Pt i - Ei so correlation 
is intriguing. As in the case of short GRBs, these evidences 
show the potential use of the E Pi i - Ei so plane to distin- 
guish among different sub-classes of GRBs, and have impor- 
tant implications for GRB/XRF/SN unification models. The 
most common explanations assume that the peculiarity of 
these events is due to pa rticular and uncommon viewing an- 
gles, a s proposed e.g. by ,YamazakL^onetokii& i Nakamura 
12003ft for GRB980425 and iRamirez-Ruiz et alJ J2005I) for 
GRB031203. Based on relativistic beaming and Doppler ef- 
fects and the assumption of a uniform jet, they find that E Pi i 
oc 5 and -Eiso oc <5 3 , where 8 is the relativistic Doppler factor 
(see Section 5.2). For large off-axis viewing angles the differ- 
ent dependence of E Pt i and Ei so on 5 would cause significant 
deviations form the E Pt i - Ei so cor relation. An a l ternat ive 
explanation has been suggested b y |Pado_fc_TJa3j j2005j) i n 
the framework of the CB model ("P ar fc De Ruiulall2004) . 
In this scenario, the vV v spectra of GRBs are characterized 
by two peaks, one at sub-MeV energies, the normal peak fol- 
lowing the E P] i - Ei ao correlation, and one in the GeV-TeV 
range. When a GRB is seen very off-axis, the same rela- 
tivistic Doppler and beaming effects discussed above would 
shift the high energy peak at low energies, making it to 
be confused with the normal low energy GRB peak. Very 
recently, off-axis scenarios have been seriously challanged 
by the detection of a very close (2 = 0.0331) and sub- 
energetic event, GRB060218/SN2006aj, which, contrary to 
GRB980425, is characterized by E P] i and _Ej BO values fully 
consistent with the E p ,j - E lao correlation dAmati et alJ 
l200dlGhisellini et al.ll2006Tl. 

Finally. iBosniak et al.1 J2006bl) . based on BATSE data of 



GRBs with unknown redshift found evidence that a part 
of GRB with indications of association with a SN are in- 
consistent with the E p ,j - Ej ao correlation, as a l so the lag- 
luminosity relation dNorris. Marani fc BonnelH |2000F) . for 
any value of z, which would confirm the existence of a pecu- 
liar class of sub-energetic, SN-associated events. See, how- 
ever, next Section for a discussion of this method. 



6 THE E p ,i - E iso CORRELATION AND GRBS 
WITH UNKNOWN REDSHIFT 

6.1 Pseudo— redshifts, GRB luminosity function, 
cosmology 

The existence of correlations between intrinsic properties of 
GRBs, emerged in the last years thanks to the increasing 
number of GRBs with known redshift, has stimulated their 
use for the estimate of pseudo-redshifts for large samples 
of BATSE GRBs. In turn, the pseudo-redshifts estimates 
have been used to compute the luminosity of large samples 
of GRBs and infer their luminosity function. This has been 
done mainly with the co rrelation between spectral l ag an d 
luminosity discovered bv lNorris. Marani fc Bonnell| (2000). 
the variability - peak luminosity correlation jRei chart^ tHil] 
l200ll iReichart fc NvsewandeJ bood iGuidorzi et alJl2005ll. 
and t he E Pr i - peak luminosity correlation iYonetoku et alJ 
2004). In addition, the E P: i - E~, and E Pi \ - E\ ao - tb cor- 
relations, given their low dispersion, have been used for the 
estima t e of cosmological p arame te rs, e.g.. iGhirlanda et alJ 
J2004I): iDai. Liang fc Xul (120041) : iLiang fc Zhand J2005fi : 
IXu. Dai fc Liand 1200 5) . in a way similar to that used with 
type la SNe . With respect to all the above correlations, 
the _B p ,i - 13i so correlation is based on a much larger sam- 
ple of GRBs, as shown in this work, and is very highly 
significant. Also, differently from the E P:1 - E 7 and E Pi i 

- Ei so - tb correlation, which show a lower dispersion but 
an higher number of possible outliers (as discussed above), 
the E Pt i - Eiso correlation does not require the detection 
of a break in the afterglow light curve nor a modelization 
of the GRB jet and circum-burst environment (as needed 
for the E Pi i - E 7 correlation). Thus, in principle the E Pi i 

- Ei so can provide the most reliable pseudo-redshift esti- 
mates. The most straightforward method is to take the flu- 
ence and spectral parameters of a GRB and compute E p ,j 
and E \so , following the same methods used bv | Amati et alJ 
J2002h or IGhirlanda. Ghisellini fc Lazzatil J2004I) . for a grid 
of 2 values (say from 0.01 to 50). The pseudo-redshift range 
will then be given by the values of z for which, accounting 
for the uncertainties on E Pt i and Ei so and for the extra- 
Poissonian dispersion of the correlation (as quantified in 
Section 4 and discussed above), the corresponding (E Pi i 
,Eiso ) points are consistent with the E p ^ - Ei so correla- 
tion at a given level of significance. In practice, if K and 
m are the normalization and index of the power-law best 
fitting the E Pi i - Ei so correlation (~100 and ~0.5 if we as- 
sume the values determined in Section 4 by accounting for 
sample variance), E p ,i and Ej so are the intrinsic peak en- 
ergy and the isotropic radiated energy at the redshift z, 
the significance of the deviation of the (E P ,i ,-Eiso ) point 
from the correlation is given by A/ yj o\ + <J%orr, where 
A = log(E P: i)~log(K)—mxlog(Ei SO ), <ta is the uncertainty 
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on A computed from oi og E 4 and oi og E iBO by error propa- 
gation and Ucorr is the extra-Poissonian dispersion of the 
correlation (based on the sample variance analysis reported 
in Section 4, one can for instance assume <r corr =0.15). Off 
course, to be more accurate one should also take into ac- 
count the uncertainties on K, m and a CO rr', however, K and 
m are correlated and have low uncertainties (Table 3) and 
for Ocorr one can conservatively assume the upper bound of 
the 90% c.l. interval (0.17). This method, as can be seen 
by testing it on GRBs with known z and published spectral 
parameter s, like e.g. th e -Bepp oSAX events included in the 



sample of lAmati et al. J2002J) or the HETE- 2 events ana- 



lyzed bv lBarraud et alJ J2003) and lSakamoto et alj l l2005rjft . 
provides reliable but often large ranges of pseudo-z (or only 
lower limits) . More precise z estimates can be obtained with 
redshift indicators partially b ased on the E p .j - Ei so correla- 
tion, like the one developed by Attcia (2003), which provides 
redshift estimates accurate to a factor of ~2 and is currently 
used to estimate pseudo-redshifts of HETE-2 GRBs, or its 
refined version proposed very recently by IPelangeon et alJ 
(2006). A caution on the use of these redshift estimators 
comes from the fact that they are partially empirical and 
thus are not supported by a complete understanding of the 
underlying physics. 

6.2 Tests and selection effects 

As already pointed out by lAmati et al] <l2002l) . given the 
relevance of the E P] i - Ei so correlation, attention has to be 
paid to the possible impact of selection effects. Recently, 
two r esearch groups jNakar fc Piranll2005l iBand fc Preecel 
120051) . by analyzing BATSE GRBs without known redshift, 
inferred that ~half jNakar fc Piranl 12005ft or even ~80% 
jBand fe Preecd 12005 ) of the whole GRB population can- 
not satisfy the correlation for any values of redshift. Thus, 
they conclude that strong selection effects are introduced in 
the various steps leading from GRB detection to the final 
z estimate and that we are measuring the redshift of only 
those events that follow the correlation. However, there are 
increasing evidences that the possible impact of selection 
effects on the E P] i - Ei ao correlation and the possible num- 
ber of outliers are much lower than argued by these au- 
thors. First, their co nclusions have been questioned by sev- 
eral other authors |G hirlanda 1 Ghiselli ni fc Firmanil 120051: 
iBosniak et al]Eo06al : IPizzichhii^tHdTlljOOfil) . who found in- 
stead that the E Pi0 bs and fluence values of most BATSE 
GRBs with unknown redshift are fully consistent with the 
E P] i - Ei ao correlation. The main source of discrepancy be- 
tween these two different conclusions lies in accounting or 
not for the observed dispersion of the correlation and for 
the uncertainties in the -Ep, bs and fluence values. When ac- 
counting for both in the way discussed in previous Section, 
it can be found that only a very small fraction of BASTE 
GRBs with unknown reds hift may be considered outliers to 
the E u .i - Ej ao correlation Ghirlanda. Ghisellini fc Firmanil 



120051: fPizzichini et al] Eo06). It has also to be noticed that 
some of these authors assumed pseudo-redshifts derived 
by other correlations, like e.g. the lag-luminosity relation, 
which are necessarily not the z values corresponding to the 
(E p ,i ,Ei a o ) with the minor deviation from the E p ,i - Ei so 
relation found with GRBs with known redshift. An impor- 
tant effect that should be taken into account when using 



BATSE data to test the correlation is that, as discussed 
in previous Section, given to the lack of coverage of the 
X-ray band below ~25 keV, BATSE is likely to overesti- 
mate t he E p ,j values at least for a fract i on of GRBs. And 
indeed, iGhirlanda. Ghisellini fc Firmanil (120051) . by fitting 
the (£ p ,i ,E ia o ) points of 442 BATSE GRBs with pseudo- 
redshifts derived by using the lag-luminosity correlation, 
find a slope and dispersion consistent with the one obtained 
with GRBs with known redshift, but a higher normaliza- 
tion, which could indicate a systematic overestimate of E p .i 
by BATSE with respect to BeppoSAX and HETE-2. I also 
note that the pos sible existence of an E u .j - Ej s o corr elation 
was suggested bv lLlovd. Petrosian fc Mallozzil i2000l) based 
on BATSE data of GRBs without known redshift. 

Secondly, as a check, I applied the pseudo-redshift es- 
timate method described above to a sample of 46 HETE-2 
GRB s with sp e ctral para meters and fluences publ ished by 
iBarraud et"afl l)2003h and ISakamoto et alj (|2005bh . I con- 



sidered only those events with both a and -E p , bs con- 
strain ed; I also disch a rged t he few events in the sam- 
ple of ISakamoto et al.1 l|2005bfl with j3>-2. When the re- 
ported best fit model is a cut-off power-law, I assumed 
for (3 the typical value of —2.5 . For those events con- 
tained in both samples. I took the spectral parameters from 
ISakamoto et alJ i2005bl) . given that it reports results based 
on dat a from both WXM and FREGATE, whereas the anal- 
ysis of IBarraud et al.1 l|2003fl is based on FREGATE data 
only. The errors on E P: i and -Ei so where derived basing on 
the errors on spectral parameters and fluences reported in 
the two references. In the estimate of _Ei so as a function of 
z I assumed a standard cosmology with fl m — 0.3, Oa=0.7 
and Ho=65 km s _1 Mpc -1 ; for the E Pi i - Ei BO correlation 
I assumed K=100, m=0.5 and ir corr =0.15 . I find that 38 
events (~83%) are consistent within la with the correlation 
(i.e. show A/y/<T^ + Ocorr ^1), 5 events (~11%) are consis- 
tent within 2a and 2 (~4%) within 3a (with values of 2.1 and 
2.7). The only event with a substantial deviation from the 
E p ,i - E iBO correlation (3.82cr) is the short GRB020531. It is 
also important to note that the pseudo-redshift ranges ob- 
tained are fully consistent with the observed z distribution. 
These results show that at least most of HETE-2 GRBs with 
unknown z are potentially fully consistent with the E P: i - 
Ei S o correlation holding for long GRBs and give a further ev- 
idence that short GRBs do not follow this correlation. Thus, 
if there are selection effects in the sample of HETE-2 GRBs 
with known redshift, on which the E Pi i - Ei so correlation is 
partly based, they are more likely due to detectors sensitiv- 
ity as a function of energy than to the subsequent processes 
leading to redshift estimate. However, the fact that the dis- 
tribution of these GRBs in the fluence-E Pi0 b s plane is con- 
sistent with that of BATSE events (Figure 16 of Sakamoto 
et al. 2005b) indicates that the possible inconsistency of a 
fraction of BATSE GRBs with unknown redshift with the 
_E Pi i - Ei S o correlation may be due to an overestimate of 
E P: obs , as a consequence of the effects discussed above and 
in Section 5.3 . 

A third important issue concerns the fact that, as discussed 
in Section 4 and shown in Figures 2 and 3, all Swift GRBs 
with known redshift are consistent with the E Pt i — Ei ao corre- 
lation. This is a strong evidence against the existence of rel- 
evant selection effects in the updated sample of GRBs with 
known redshift on which the E Pt i - Ei so correlation study 
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here presented is based, because: a) the burst detection sen- 
sitivity of Swift /BAT in ~ 15-300 keV is comparable to that 
of BATSE and better than that of BeppoSAX and HETE-2 
(see however Band 2003 for a comparison of the sensitivities 
of these different detectors as a function of energy), which 
in the past years contributed nearly all the measurements 
on which the E Pt i - Ei so correlation was based; b) the very 
fast and precise localization capabilities of Swift /XRT al- 
lowed to substantially reduce the selection effects also in the 
process leading to redshift estimate (GRB precise localiza- 
tion, optical follow-up, optical afterglow and/or host galaxy 
detection and spectroscopy). A drawback of BAT is that it 
can provide an estimate of -E Pl obs only for a small fraction 
of events (15-20% at most). 

Basing on the above, it is unlikely that the E Pt i - Ei so corre- 
lation is strongly affected by selection effects. Anyway, the 
existence of sub-classes of GRBs not following it cannot 
be excluded. As discussed in previous Section, the analy- 
sis of the sample of GRBs with known redshift and i?p,obs 
shows that, in addition to short GRBs, a class of sub- 
energetic events (like GRB980425 and possibly GRB031203) 
with spectral-energy properties inconsistent with the corre- 
lation may exist. The possibile existence of a fraction of 
long GRBs is also predic ted by some GR B synthesis mod- 
els, like, e.g., the one by iDonaghvl j200d) based on off-jet 
relativistic kinematics effects. Obviously, the most reliable 
test of the E P] i - i?j so correlation and of the existence of 
one or more sub-classes of outliers will come from the on- 
going quick enlargement of the sample of GRBs with known 
redshift and i? p , bs both in number, thanks to Swift fast 
and precise localizations, and in the coverage of the -E p>0 b s 
-fluence plane, as allowed by the GRB experiments with dif- 
ferent energy bands and sensitivity presently operating, like 
Swift /BAT, HETE-2, Konus/Wind, Suzaku/WAM, INTE- 
GRAL/ISGRI, or that will fly in the next future, like, e.g., 
those on board AGILE and GLAST. 
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